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Abstract 

In this paper, the lowest order constrained variational (LOCV) method has been used to inves- 
tigate the magnetic properties of spin polarized neutron matter in the presence of strong magnetic 
field at zero temperature employing AVis potential. Our results indicate that a ferromagnetic 
phase transition is induced by a strong magnetic field with strength greater than 10 18 G, leading 
to a partial spin polarization of the neutron matter. It is also shown that the equation of state 
of neutron matter in the presence of magnetic field is stiffer than the case in absence of magnetic 
field. 
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I. INTRODUCTION 



The magnetic field of neutron stars most probably originates from the compression of 
magnetic flux inherited from the progenitor star Using this point of view, Woltjer has 
predicted a magnetic field strength of order 10 15 G for the neutron stars |2j. The field can 
be distorted or amplified by some mixture of convection, differential rotation and magnetic 
instabilities {3, 4|. The relative importance of these ingredients depend on the initial field 
strength and rotation rate of the star. For both convection and differential rotation, the 
field and its supporting currents are not likely to be confined to the solid crust of the star, 
but distributed in most of the stellar interior which is mostly a fluid mixture of neutrons, 
protons, electrons, and other more exotic particles Thompson et al. |5| argued that 
the newborn neutron stars probably combine vigorous convection and differential rotation 
making it likely that a dynamo process might operate in them. They expected fields up 
to 10 15 — 10 16 G in neutron stars with few millisecond initial periods. On the other hand, 
according to the scalar virial theorem which is based on Newtonian gravity, the magnetic 
field strength is allowed by values up to 10 18 G in the interior of a magnetar However, 
general relativity predicts the allowed maximum value of neutron star magnetic field to be 
about 10 18 — 10 20 G [7]. By comparing with the observational data, Yuan et al. [8] obtained 
a magnetic field strength of order 10 19 G for the neutron stars. 

The strong magnetic field could have important influences on the interior matter of a 
neutron star. Many works have dealt with study of the magnetic properties and equation 



of state of the neutron star matter M 



20j and quark star matter 2ll426| in the presence of 



strong magnetic fields. Some authors have considered the influence of strong magnetic fields 



on the neutron star matter within the mean field approximation 



Yuan et al. 



9] using 



the nonlinear a — u> model, showed that the equation of state of neutron star matter becomes 
softer as the magnetic field increases. Also, Broderick et al. [hJ employing a field theoretical 
approach in which the baryons interact via the exchange of a — u — p mesons, observed that 
the softening of the equation of state caused by Landau quantization is overwhelmed by 
stiffening due to the incorporation of the anomalous magnetic moments of the nucleons. It 
has been shown that the strong magnetic field shifts /9-equilibrium and increases the proton 



fraction in the neutron star matter 



10-121. Yue et al. 



13] have studied the neutron star 



matter in the presence of strong magnetic field using the quark-meson coupling (QMC) 
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model. Their results indicate that the Landau quantization of charged particles causes a 
softening in the equation of state, whereas the inclusion of nucleon anomalous magnetic 
moments lead to a stiffer equation of state. The effects of the magnetic field on the neutron 

n 

star structure, through its influence on the metric has been studied by Cardall et al. |27J. 
Their results show that the maximum mass, in a static configuration for neutron star with 
magnetic field, is larger than the maximum mass obtained by uniform rotation. Through 
a field theoretical approach (at the mean field level) in which the baryons interact via the 
exchange of a — u — p mesons, Broderick et al. [ijj have considered the effects of magnetic 
field on the equation of state of dense baryonic matter in which hyperons are present. They 
found that when the hyperons appear, the pressure becomes smaller than the case of pure 
nucleonic matter for all fields. Within a relativistic Hartree approach in the linear a — u; — p 
model, the effects of magnetic field on cold symmetric nuclear matter and the nuclear matter 
in /3-equilibrium have been investigated by Chakrabarty et al. 15]. Their results suggest 
that the neutron star mass is practically insensitive to the effects of the magnetic fields, 
whereas the radius decreases in intense fields . 

In some studies, the neutron star matter was approximated by a pure neutron matter. 
Isayevet a,. Q conned the neutron ^ter in a strong magnetic field with the Skyrme 
effective interaction and analyzed the resultant self-consistent equations. They found that 
the thermodynamically stable branch extends from the very low densities to the high den- 
sity region where the spin polarization parameter is saturated, and neutrons become totally 



spin polarized. Perez-Garcia et al. |18M20| studied the effects of a strong magnetic field 
on the pure neutron matter with effective nuclear forces within the framework of the non- 
relativistic Hartree-Fock approximation. They showed that in the Skyrme model there is 
a ferromagnetic phase transition at p ~ 4po(po = 0.16/m -3 is the nuclear saturation den- 
sity), whereas it is forbidden in the DIP model Beside these, they found that the 
neutrino opacity of magnetized matter decreases compared to the nonmagnetized case for 
the magnetic field greater than 10 17 G 19]. However, more realistically, for the problem of 
the neutron star matter in astrophysics context, it is necessary to consider the finite tern- 

^ n nn n 

24j and finite proton fraction effects |9_Hl5j. Isayev et al. [17] have 



perature 
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shown that the influence of finite temperatures on spin polarization remains moderate in 
the Skyrme model, at least up to temperatures relevant for protoneutron stars. It has been 
also shown that for SLyA effective interaction, even small admixture of protons to neutron 
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matter leads to a considerable shift of the critical density of the spin instability to lower 
values. For SkI5 force, however, a small admixture of protons to neutron matter does not 
considerably change the critical density of the spin instability and increases its value 28 1. 
In our previous works, we have studied the spin polarized neutron matter |29| , symmetric 



nuclear matter [30], asymmetric nuclear matter [31], and neutron star matter [3l| at zero 
temperature using LOCV method with the realistic strong interaction in the absence of 
magnetic field. We have also investigated the thermodynamic properties of the spin polarized 



neutron matter |32J, symmetric nuclear matter [33J, and asymmetric nuclear matter [34| at 
finite temperature with no magnetic field. In the above calculations, our results do not show 
any spontaneous ferromagnetic phase transition for these systems. In the present work, we 
study the magnetic properties of spin polarized neutron matter at zero temperature in the 
presence of the strong magnetic field using LOCV technique employing AV\s potential. 

II. LOCV FORMALISM FOR SPIN POLARIZED NEUTRON MATTER 

We consider a pure homogeneous spin polarized neutron matter composed of the spin-up 
(+) and spin-down (— ) neutrons. We denote the number densities of spin-up and spin-down 
neutrons by p^ and p^~\ respectively. We introduce the spin polarization parameter (5) 
by 

p(+) _ o(-) 

where — 1 < 8 < 1, and p = p^ + > + p*"- 1 is the total density of system. 

In order to calculate the energy of this system, we use LOCV method as follows: we 
consider a trial many-body wave function of the form 

V> = F<j>, (2) 

where is the uncorrelated ground-state wave function of N independent neutrons, and F is 
a proper iV-body correlation function. Using Jastrow approximation |35|, F can be replaced 
by 

F = SUf(ij), (3) 

i>j 
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where S is a symmetrizing operator. We consider a cluster expansion of the energy functional 
up to the two-body term, 

rnm 1 { ^ m) F I F U) 
Em -N^4W- El + ET (4) 

Now, we calculate the energy per particle up to the two-body term for two cases in the 
absence and presence of the magnetic field in two separate sections. 



A. Energy calculation for the spin polarized neutron matter in the absence of 
magnetic field 

The one-body term E\ for spin polarized neutron matter in the absence of magnetic field 
(B = 0) is given by 

(B=0) = y Shk^p^ 

where kp = (Qtc 2 p^)^ is the Fermi momentum of a neutron with spin projection i. 
The two-body energy E2 is 

1 



where 



Ei 2 B=0) = ^E^>( 12 )IU-J2>, (6) 



"(12) = -|^[/( 12 )> ^ /(I 2 )]] + /(12)V(12)/(12). 



In the above equation, /(12) and 1^(12) are the two-body correlation function and nuclear 



potential, respectively. In our calculations, we employ the AVi$ two-body potential 



36|, 



18 

l(p) 



V(12) = E v(p) (^)°il (7) 
P =i 



where 

(p_i-i8) = i )(7l cr2)Tl T2) (cri.o- 2 )(ri.r 2 ), S"^, 5 , i 2 (ri.r 2 



1 

L.S, L.S(ti.t 2 ), L 2 , L 2 ((Ti.o- 2 ), L 2 (ri.r 2 ), L 2 (cri.o- 2 )(ri.r 2 ), 
(L.S) 2 , (L.S) 2 (ri.r 2 ), T 12 , (a 1 . ( r 2 )T 12 , S 12 T 12 , (r 2l + r z2 ). (8) 



In the above equation, 

S12 = [3(cri.f)(o- 2 .f) - cri.cr 2 ] 



is the tensor operator and 

T12 = [3(ri.r)(r 2 .f)-ri.r 2 ] 

is the isotensor operator. The above 18 components of the AVi% two-body potential are 
denoted by the labels c, a, r, err, t, tr, Is, 1st, 12, 12a, I2r, I2ar, ls2, Is2t, T, oT , tT, and 



rz, respectively 



36l | . In the LOCV formalism, the two-body correlation function /(12) is 



considered as follows 



37|, 



3 

1 

k=l 



/(12) = E/^MPf, 



where 



P 



(fc=l-3) 



: ( 2 ) 



11 

12 - l 4 4 - iz v 2 ■ g- iz ■ g 



1 :f' 

6 



(9) 



(10) 



The operators O l2 and 0} 2 are given in Eq. (jSJ). Using the above two-body correlation 
function and potential, after doing some algebra, we find the following equation for the 
two-body energy: 



E (B=0) 



11 * 2 

— CTzl - cr z o SS Z 
2 2 1 



A fill) V p _ (_!)«*,.] 

/f T 0,771 9 

<fr[ (V) +-^-({K-3^ + K-3K T + 2(y T -3y CTT )-2y r ja« (fy 

+ [Vl 2 -3y /2(T + ^ 2T -3^ 2(Tr ]ci 1 ) 2 (V))(^ 1) ) 2 l+ E U fc) ' 2 «i fc)2 (V) 



r 



k=2,3 



+^r({Vc + V a + V T + V aT + (-6k + U)(V tT + V t ) — (k— l)(V lST + V u ) 
+2[V T + V aT + (-6* + U)Vrr - V TZ }}a^\k } r) 
+ [V l2 + Vna + V l2T + V l2aT ]c^\k f r) + [V ls2 + V ls2T }d^\k f r))f^ 
+f?Ok + V lsr - 2(V l2 + V l2a + V l2uT + V l2T ) - 3(V ls2 + V ls2T )}b 2 a {k f r)fWf^ 



(11) 



where a = {J, L, S, S z } and the coefficient a^ 1 -* 2 , etc., are defined as 



(12) 



a a ] ( X ) = 3?\PIj-l,S,( x ) + 7-^+1,5,0*0], 



aj> (x) = x 2 [ 7 /j-i,5 z (x) + PI J+ll s.(x)], 



(13) 
(14) 
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with 



and 



b^\x) = x 2 [f3 23 Ij^ Sz (x) - (3 23 I J+1 , Sz (x)}, 

c£ )2 (x) = x 2 vJ L>Sz {x), 
c£ )2 (x) = x 2 [r} 2 Ij- hSz (x) + v 2 I J+ i,s z (x)\, 

cf 2 {x) = x 2 [r] 3 Ij_ hSz (x) + v 3 Ij+i,s z ( x )]i 
d^ 2 (x) = x 2 [^Ij^, Sz (x) + X 2 Ij +1 , Sz (x)i 
Sf(x) = x 2 [i 3 Ij_ 1:Sz {x) + X 3 Ij +1 , Sz (x)}, 



' 2J+1' 7 2J+1' P23 2J+1 ' 

TIT ^ J 2 {J + l) J 3 + 2J 2 + 3J + 2 
* = L(L + 1), ^ = ^7^, -3 = ^ , 

J(J 2 + 2J+1) J{J 2 + J + 2) 

V2 ~ 2jTl ' V3 ~ 2J+1 ' 

_ J 3 + 2 J 2 + 2 J + 1 _ J(J 2 + J + 4) 
^ 2 ~ 2J + 1 ' ^ 3 ~ 2J+ 1 ' 

J(J 2 + J+1) _J 3 + 2J 2 + 5J + 4 
A2 " 2J+1 ' As " 27TI ' 



Ij,S z (x)= Jdq q 2 PsMJ 2 {xq)- 



In the last equation Jj(x) is the Bessel function and Ps z (q) is defined as 
PsM) = ^[(^ 1 ) 3 + (^ 2 ) 3 -^((^ 1 ) 2 + (^ 2 ) 2 )g 

- ^m 1 ) 2 - (k^m- 1 + q 3 } 

for - k a / 2 \ < q < Hk' 7 / 1 + k ff F "\, 

Ps z {q) = \™in{{k^)\{krf) 

for q < \\k F zl - fcj,* a |, and 

PsM = 

for g > ||/c^ zl + kp z2 \, where o z \ or a z2 = +1, —1 for spin up and down, respectively. 



B. Energy calculation of spin polarized neutron matter in the presence of magnetic 
field 



Now we consider the case in which the spin polarized neutron matter is under the influence 
of a strong magnetic field. Taking the uniform magnetic field along the z direction, B = Bk, 
the spin up and down particles correspond to parallel and antiparallel spins with respect to 
the magnetic field. Therefore, the contribution of magnetic energy of the neutron matter is 

E M = —M Z B, (30) 

where M z is the magnetization of the neutron matter which is given by 

M z = NfiJ. (31) 

In the above equation, //„ = —1.9130427(5) is the neutron magnetic moment (in units of 
the nuclear magneton). Consequently, the energy per particle up to the two-body term in 
the presence of magnetic field can be written as 

E([f}) = E{ B=0) + E { 2 B=0) - ^ n B5, (32) 

where e[ B=0) and =0) are given by Eqs. (jSJ) and (1111) . respectively. It should be noted 
that in usual thermodynamic treatments the external magnetic field energy (g- J_dV B 2 ) is 
usually left out since it does not affect the thermodynamic properties of matter [38(. In fact 
the magnetic field energy arises only from the magnetostatic energy in the absence of matter, 
but we are interested in the contribution of internal energy which excludes the energy of 
magnetic field. Therefore, the magnetic field contribution, E mag = |^-, which is the energy 
density (or "magnetic pressure") of the magnetic field in the absence of matter is usually 



omitted 
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Now, we minimize the two-body ener gy w ith respect to the variations in the function /W 
subject to the normalization constraint 



N 



£<v>l,-/ 2 (i2) y)- = o 



(33) 



where in the case of spin polarized neutron matter, the function hs z {r) is defined as follows, 



hsA r ) 



1-9 



21 "1/2 



S z = ±l 



(34) 



S z = 



From minimization of the two-body cluster energy, we get a set of coupled and uncoupled 



differential equations which are the same as those presented in Ref . [39| , with the coefficients 
replaced by those indicated in Eqs. (TT2"j) — ( 12"U|) . By solving these differential equations, we 
can obtain correlation functions to compute the two-body energy. 



III. RESULTS AND DISCUSSION 



Our results for the energy per particle of spin polarized neutron matter versus the spin 
polarization parameter for different values of the magnetic field at p = 0.2 fm~ 3 have been 
shown in Fig. [IJ We have found that for the values of magnetic field below 10 18 G, the 
corresponding energies of different magnetic fields are nearly identical. This shows that 
the effect of magnetic field below B ~ fO 18 G is nearly insignificant. From Fig. [H we 
can see that the spin polarization symmetry is broken when the magnetic field is present 
and a minimum appears at — f < 5 < 0. By increasing the magnetic field strength from 
B ~ fO 18 G to B ~ fO 19 G, the value of spin polarization corresponding to the minimum 
point approaches — f. We also see that by increasing the magnetic field, the energy per 
particle at minimum point (ground state energy) decreases, leading to a more stable system. 
For each density, we have found that above a certain value of the magnetic field, the system 
reaches a saturation point and the minimum energy occurs at 5 = —1. For example at 
p = 0.2 fm~ 3 , for B > f .8 x fO 19 G, the minimum energy occurs at 5 = — 1 . However, this 
threshold value of the magnetic field increases by increasing the density. In Fig. [2j we have 
presented the ground state energy per particle of spin polarized neutron matter as a function 
of the density for different values of the magnetic field. For each value of the magnetic field, 
it is shown that the energy per particle increases monotonically by increasing the density. 
However, the increasing rate of energy versus density increases by increasing the magnetic 
field. This indicates that at higher magnetic fields, the increasing rate of the contribution 
of magnetic energy versus density is more than that at lower magnetic fields. In order to 
clarify this behavior, we have presented the energy contribution of spin polarized neutron 
matter up to the two-body term in the cluster expansion (E\ + E2), and the magnetic energy 
contribution (Em) separately, as a function of density in Fig. |3j This figure shows that for 
the spin polarized neutron matter, the difference between the magnetic energy contributions 
(Em) of different magnetic fields is substantially larger than that for the energy contribution 
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(E1 + E2). Fig. H] shows the ground state energy per particle of spin polarized neutron matter 
as a function the magnetic field for different values of density. We can see that by increasing 
the magnetic field up to a value about 10 18 G, the energy per particle slowly decreases, 
and then it rapidly decreases for the magnetic fields greater than this value. This indicates 
that above B ~ 10 18 G, the effect of magnetic field on the energy construction of the spin 
polarized neutron matter becomes more important. 

In Fig. [51 the spin polarization parameter corresponding to the equilibrium state of the 
system is plotted as a function of density for different values of the magnetic field. It is 
seen that at each magnetic field, the magnitude of spin polarization parameter decreases 
by increasing the density. Fig. also shows that for the magnetic fields below 10 18 G, at 
high densities, the system nearly becomes unpolarized. However, for higher magnetic fields, 
the system has a substantial spin polarization, even at high densities. In Fig. |6l we have 
plotted the spin polarization parameter at the equilibrium as a function of the magnetic 
field at different values of density. This figure shows that below B ~ 10 18 G, no anomaly is 
observed and the neutron matter can only be partially polarized. This partial polarization 
is maximized at lower densities and amounts to about 14% of its maximum possible value of 
— 1 . From Fig. [61 we can also see that below B ~ 10 17 G, the spin polarization parameter 
is nearly zero. This clearly confirms the absence of the magnetic ordering for the neutron 
matter up to B ~ 10 17 G. For the magnetic fields greater than about 10 18 G, it is shown that 
the magnitude of spin polarization rapidly increases by increasing the magnetic field. This 
shows a ferromagnetic phase transition in the presence of a strong magnetic field. For each 
density, we can see that at high magnetic fields, the value of spin polarization parameter 
is close to —1. The corresponding value of the magnetic field increases by increasing the 
density. 

The magnetic susceptibility (x) which characterizes the response of a system to the 
magnetic field, is defined by 



In Fig. [7J we have plotted the ratio x/iV|/i n | for the spin polarized neutron matter versus 
the magnetic field at three different values of the density. As can be seen from Fig. [7J for 
each density, this ratio shows a maximum at a specific magnetic field. This result confirms 
the existence of the ferromagnetic phase transition induced by the magnetic field. We see 




(35) 
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that the magnetic field at phase transition point, B m , depends on the density of the system. 
Fig. |8] shows the phase diagram for the spin polarized neutron matter. We can see that by 
increasing the density, B m grows monotonically. It explicitly means that at higher densities, 
the phase transition occurs at higher values of the magnetic field. 

From the energy of spin polarized neutron matter, at each magnetic field, we can evaluate 
the corresponding pressure (Pkinetic) using the following relation, 



P k >neUc(p,B)=p*(^^l) (36) 



{ dp 



Our results for the kinetic pressure of spin polarized neutron matter versus the density for 
different values of the magnetic field have been shown in Fig. [9j It is obvious that with 
increasing the density, the difference between the pressure of spin polarized neutron matter 
at different magnetic field becomes more appreciable. Fig. [9] shows that the equation of 
state of the spin polarized neutron matter becomes stiffer as the magnetic field strength 
increases. This stiffening is due to the inclusion of neutron anomalous magnetic moments. 
Th,s i, in agreement with the results obta,„ed in Refs. QQ. It 8 hou.d be noted he r e 
that to find the total pressure related for the neutron star structure, the contribution from 
the magnetic field, P m a g = f^-, should be added to the kinetic pressure Q, [l4|. However, 
in this work we are not interested in the neutron star structure and have thus omitted the 
contribution of "magnetic pressure" in our calculations for neutron matter [16] . This term, 
if included, simply adds a constant amount to the curves depicted in Fig. [9j 



IV. SUMMARY AND CONCLUDING REMARKS 

We have recently calculated several properties of the spin polarized neutron matter in 
the absence of magnetic field using the lowest order constrained variational method with 
AV±s potential. In this work, we have generalized our calculations for spin polarized neutron 
matter in the presence of strong magnetic field at zero temperature using this method. We 
have found that the effect of magnetic fields below B ~ 10 18 G is almost negligible. It 
was shown that in the presence of magnetic field, the spin polarization symmetry is broken 
and the energy per particle shows a minimum at — 1 < 5 < 0, depending on the strength 
of the magnetic field. We have shown that the ground state energy per particle decreases 
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by increasing the magnetic field. This leads to a more stable system. It is seen that the 
increasing rate of energy versus density increases by increasing the magnetic field. Our 
calculations show that above B ~ 10 18 G, the effect of magnetic field on the properties of 
neutron matter becomes more important. In the study of spin polarization parameter, we 
have shown that for a fixed magnetic field, the magnitude of spin polarization parameter at 
the minimum point of energy decreases with increasing density. At strong magnetic fields 
with strengths greater than 10 18 G, our results show that a field-induced ferromagnetic 
phase transition occurs for the neutron matter. By investigating the magnetic susceptibility 
of the spin polarized neutron matter, it is clear that as the density increases, the phase 
transition occurs at higher values of the magnetic field. Through the calculation of pressure 
as a function of density at different values of the magnetic field, we observed the stiffening 
of the equation of state in the presence of the magnetic field. 

Finally, we would like to address the question of thermodynamic stability of such neutron 
stars at ultra-high magnetic fields. One may wonder if the effect of magnetic pressure, 
Pmag = f^-, which we have omitted here, is added to the kinetic pressure Pkinetic, then 
at ultra-strong magnetic fields, the system might become gravitationally unstable due to 
excessive outward pressure. For the fields considered in this work (up to 10 20 G), this 
scenario does not seem likely j7|. We note that the increase of magnetic field leads to 
stiffening of the equation of state (Fig. [9]) which in turn leads to larger mass and radius for 



the neutron star 40j. This in turn increases the effect of gravitational energy, offsetting the 
increased pressure. We also note that the existence of a well-defined thermodynamic energy 
minimum for all fields considered in our work indicates the thermodynamic stability of our 
system. The existence of such well-defined minimum energy is unaffected by the addition of 
magnetic energy. The detailed analysis of such situations along with accompanying change 
in proton fraction is a possible avenue for future research. 
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FIG. 1: The energy per particle versus the spin polarization parameter (5) for different values of 
the magnetic field (B) at p = 0.2 /m -3 . 
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FIG. 2: The ground state energy per particle as a function of the density at different values of the 
magnetic field (B). 
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FIG. 3: The energy contribution of spin polarized neutron matter in the cluster expansion up to the 
two body term (E\ + E2) for the magnetic fields B = 10 18 G (solid curve) and B = 10 19 G (dashed 
dotted curve), and the contribution of magnetic energy (Em) for magnetic fields B = 10 18 G 
(dashed curve) and B = 10 19 G (dashed dotted dotted curve). 
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FIG. 4: The ground state energy per particle as a function of the magnetic field (B) at different 
values of the density (p) . 
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FIG. 5: The spin polarization parameter at the equilibrium state of the system as a function of 
the density at different values of the magnetic field (B) . 
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FIG. 6: The spin polarization parameter corresponding to the equilibrium state of the system as 
a function of the magnetic field (B) at different values of the density (p). 
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FIG. 7: The magnetic susceptibility (x/N\(i n \) as a function of the magnetic field (B) at different 
values of the density (p) . 
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FIG. 8: Phase diagram for the spin polarized neutron matter in the presence of strong magnetic 
field. 
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FIG. 9: The equation of state of spin polarized neutron matter for different values of the magnetic 
field (B). 
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